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ABSTRACT 

Using a flexible galactic model with variable stellar velocity anisotropy, I apply the classical Jeans 
mass modeling approach to the five dwarf spheroidal galaxies with the largest homogeneous datasets 
of stellar line-of-sight velocities (between 330 and 2500 stars per galaxy) - Carina, Fornax, Leo I, 
Sculptor, and Sextans. I carry out an exhaustive model parameter search, assigning absolute proba¬ 
bilities to each parameter combination. My main finding is that there is a well defined radius (unique 
for each galaxy) where the Jeans analysis constraints on the enclosed mass are tightest, and are much 
better than the constraints at previously suggested radii (e.g. 300 pc). For Carina, Fornax, Leo I, 
Sculptor, and Sextans the enclosed DM mass is 0.94±0.20 (at 410 pc), 7.1±0.9 (at 925 pc), 1.75±0.20 
(at 390 pc), 2.59 ±0.42 (at 435 pc), and 2.3 ±0.9 (at 1035 pc), respectively (two-sigma uncertainties; 
in 10^ Mq units). Local DM density has the tightest constraints at smaller (and also unique for 
each galaxy) radii. The largest central DM density constraint is for Sculptor: po ± 0.09 Mq pc“^ 

(at two-sigma level). I show that the DM density logarithmic slope is totally unconstrained by the 
Jeans analysis at all the radii probed by the data (and not just at the center, as was demonstrated 
before). Stellar velocity anisotropy has only very weak constraints. In particular, pure central tangen¬ 
tial anisotropy is ruled out at better than two sigma level for three dwarfs, and the data are consistent 
with the global stellar velocity isotropy for all the five galaxies. 

Subject headings: galaxies: dwarf — Local Group — galaxies: kinematics and dynamics — dark 
matter — methods: numerical 


1 . INTRODUCTION 

Dwarf spheroidals (dSphs) are very faint galaxies 
(some of the recently discovered dSp h satellites of Milky 
Way are the faintest known galaxies. iMartin et al.jl200'^ 
which tend to cluster around large galaxies (iMateoll 199811 . 
Among the 50 or so Local Group dSphs discovered so 
far, only two (Tucana and Cetus) are relatively isolated; 
the rest are located in the vicinity of either Milky Way 
or M31, and are believed to be their satellites. Despite 
their dull appearances (no ongoing sta r formation; little 
or no interstellar medium; lMateol[l998ll , dSphs are fasci¬ 
nating objects - largely due to the fact that they contain 
significant quantities of dark matter (DM), and as such 
represent the smallest observed scale of DM clustering. 

State of the art cosmological simulations predict that 
large galaxies should contain many smaller sub-halos 
which have not had time to be completely destroyed by 
the tidal field of the host halo. dSphs ar e believed to 
represent at least some of this substructure (IMoore et al.l 

dSphs present two major challenges to the standard 
(ACDM) cosmological model. First, despite the recent 
advances on both theoretical and observational sides, 
there still appears to be a factor of a few discrepancy 
between the observed and predicted numbers of Milky 
Way satellites - so called “missing satellites problem” 
(|Moore et al.l I1999H . Second, some indirect evidence 
(such as the “timing p aradox” for the five globular clus¬ 
ters of Fornax dSphs; iGoerdt et al.ll200ll suggests that 
DM distribution in dSphs has a flat core, which is in 
line with what is deduced f or larger dwarf and low 
surface brightness galaxies (jvan den Bosch fc SwatersI 


20011 Ide Blok fc Bosmal 120021 : iMarchesini et al.l 120021 : 
Gentile et al.l I2005D . but is at odds with the results of 


many cosmological simulations, predicting central di- 
vergent DM cusps w ith the logarithmic slope 7 ~ 1 
(|Navarro et al.l[T997f) . Different theoretical mechanisms 
suggested to rectify this “cusp-core” problem produce 
different velocity anisotropy profiles f or stars. For exam¬ 
ple. our stel lar feedback mechanism (|Mashchenko et ^ 
l2006bl l2008fl , where the DM cusp is flattened via gravi¬ 
tational heating from the interstellar gas sloshed around 
by supernovae, results in isotropic stellar cores. On 
the other hand, th e dynamical friction mechanism of 
lEl-Zant et al.l (|2001ll . where the cusp is heated by mas¬ 
sive baryonic clumps passively spiraling in towards the 
galactic center, sho uld produce notice able central tan¬ 
gential anisotropy (jTonini et al.l I2006D . Gonstraining 
stellar anisotropy profiles in spheroidal galaxies can 
hence be a valuable indirect way of identifying the correct 
t heoretical mechan ism of DM core formation. 

lAaronsonl (|1983fl was the first who deduced the pres¬ 
ence of significant amounts of DM in a dSph (Draco), 
based on the measurements of the line-of-sight veloci¬ 
ties of only three stars. Since then, the number of stars 
in dSphs with known line-of-sight velocities has grown 
dramatically: recently published homogen eous datasets 
have u p 2500 stars per galaxy (for Fornax; iWalker et al.l 
l2009c[ f. Mass modeling of dSphs has dramatically im¬ 
proved as well: froi n simple global estimates based on 
the virial theorem (lAaronsonl 1198311 . to spatial models 
with strong simplifying assumptions (such as “mass fol¬ 
lows light”; lMated[l998f) , to more recent modeling efforts 
with fewer assumptions (Walker et al.l[2009cl : iWolf et ^ 
120091:l^rigari et ffll2008 f. 
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The classical method of dSph mass modeling is to solve 
the spherical Jeans equation for a range of models with 
different DM density profiles, and then to compare the 
predicted stellar line-of-sight velocity dispersion profiles 
with the observed ones - using either or maximum 
likelihood techniques to find the best models. The spher¬ 
ical symmetry assumption is reasonably accurate given 
the spheroidal appearances of dSphs (ellipticity e is usu¬ 
ally ~ 0.3 or less), absence of disk-li ke structures , and 
negligible rotation in these systems (jMateol 11 9981 1. To 
solve the spherical Jeans equation (that is, to derive the 
stellar velocity dispersion as a function of radius), one 
has to specify two functions: total (DM -|- stars) density, 
and stellar velocity anisotropy (/3) profiles. To reduce the 
dimensionality of the problem, one usually assumes a cer¬ 
tain shape of the /3 radial profile - either a constant, or 
the Osipkov-Merritt profile (purely isotropic at the cen¬ 
ter, purely radially anisotropic in the infinity; lOsiokoAl 
[ma IMerrittlfl^ . 

There are no reasons to believe that stellar anisotropy 
is constant in dSphs. For example, simple dynami¬ 
cal models which start with an initially non-equilibrium 
stellar configuration inside the DM halo (“cold col- 
lanse”. iMashchenko fc SilTsI 12004 and “hot collapse”, 
iMashchenko fc Sillsll2005ll tend to produce the stellar ve¬ 
locity isotropy at the center and a variable degree of 
radial anisotropy in the outskirts of the relaxed stellar 
body. As already mentioned, the central parts of dwarf 
galaxies can be significantly affected by dense and violent 
baryons, which may either rand omize stellar velocities 
(IMashchenko et al.ll2006bl 1200811 or indu ce a significant 
tangential anisotropy (iTonini et al.l[200^ . On the obser¬ 
vational side, dSphs are known to be non-homogeneous 
objects, with younger, more metal-rich and kinematically 
colder star populati ons concentrated towards the center 
((Tolstov et al.ll2?iol l: it would be strange if the stellar 
velocity anisotropy would be the same across these dif¬ 
ferent populations of stars. 

Given the expectation that /3 is not constant across 
a dSph and that the important properties of dSphs 
(such as the central logarithmic slope of the DM den¬ 
sity) are strongly degenerate with respect to the un¬ 
known stellar veloc ity anisotropy (|Strigari et al.l 120081 : 
IWalker et al.ll 200 ^ . proper mass modeling of dSphs has 
to include radially variable stellar anisotropy. To the best 
of my knowledge, the only example of a full-fledged Jeans 
mass modeling of d Sphs with variable P is the work by 
iStrigari et al.l ( 20081 1 . Here I outline the main diff erences 
of my approach with that of iStrigari et al.l (|2008ll . 

1. In this paper, I carry out mass modeling of the lim¬ 
ited number (five; Carina, Fornax, Leo I, Sculptor, 
and Sextans) of Galactic dSphs with the highest 
quality, homogeneous data: large (> 300 stars) ho- 
mogeneous datasets o f stellar line-of-sight velocities 
(jWalker et al.ll200'^ and accurate stellar surface 
brightness profiles derived in a homogeneous man- 
ner in this paper. Th is is in contrast to the work of 
IStrigari et al.l (1200811 . who analyzed a much larger 
set of Local Group dSphs, with heterogeneously de¬ 
rived data and a wide range of data quality. 

2. I employ a “brute force” optimization while search¬ 
ing for the best fitting models. This approach gen¬ 


erated a wealth of statistical data, which can be 
used for testing a variety of different hypotheses 
about the distribution of stars and DM in dSphs. 
The data are available online. 

3. The “brute force” approach coupled with a flexible 
dSph model allowed me to find the dSph param¬ 
eters which the Jeans analysis can constrain very 
well. Specifically, I found that the total enclosed 
mass in a dSph can be constrained to a high ac¬ 
curacy (two-sigma uncertainty as good as ±15%, 
for Fornax) at a certain radius, which is different 
for different dwarfs. This information can be valu¬ 
able for matching the Galactic dSphs to the pre¬ 
dictions of cosmological simulations. Similarly, I 
derived tight constraints on the local DM density 
at a certain radius (unique for each galaxy), which 
can be used in the research aimed at detecting DM 
in dSphs via its annihilation signal. 

4. Unlike IStrigari et all (|2008ll . I properly account for 
the self-gravity of stars. This can be important for 
dSphs with dense stellar cores (e.g. Sculptor). 

5. I use an advanced (with the variance due to ran¬ 
dom locations of stars within one radial bin re- 
moved; see SI4.1II model fitting approach. Unlike 
IStrigari et al.l (| 200 ^ . who employed a maximum 
likelihood technique, my approach is insensitive to 
the (unknown) shape of the line-of-sight star veloc¬ 
ity probability distribution function (PDF). 

I had two main goals for this paper. First, I wanted to 
find out which (if any) dSph parameters can be meaning¬ 
fully constrained via the classical, Jeans mass modeling 
approach, if such an analysis is pushed to the extreme 
(large homogeneous observational datasets of stellar line- 
of-sight velocities; a very flexible galactic model; careful 
numerical analysis consuming a significant - 3 x 10® cpu 
hours in my case - amount of supercomputing time; a 
“brute force” optimization which ensures that no hid¬ 
den “valleys” and local and global minima in the multi¬ 
dimensional likelihood manifold are missed). The second 
goal was to generate high-quality data which can be used 
to dramatically reduce computational time (by restrict¬ 
ing the model parameter space to be explored) in future 
post-Jeans (analyzing the full shape of the stellar velocity 
PDF) mass modeling efforts - which will be the subject 
of Paper II in this series. 

As the dynamical state of the outer parts of dSphs 
is s till a matter of de bate (specifically, some authors, 
e.g. iMunoz et al.ll2006l argue that the observed outskirts 
of these galaxies are undergoing a severe tidal disrup¬ 
tion and hence are not in dynamic equilibrium, which 
wo uld invalidate the Jeans analysis for those parts; but 
see iSegall et al.l 120071 for the opposite example), in the 
present study I put the main emphasis on recovering the 
inner structure of dSphs. Unfortunately, the most inter¬ 
esting aspect of the inner structure of these galaxies - 
the “DM cusp or core” question ~ cannot be addressed 
via the Jeans analysis alone, as the central logarithmic 
slope of the DM density profile 7 is known to be degen¬ 
erate in this method. This fo llows from the properties o f 
the spherical Jeans equation ((Binnev fc TremainallQSTll . 
and is clearly visible even in the mass modeling with a 
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constant stellar anisotropy (iWalker et alJl^OQch . In this 
paper I obtain a more general result that 7 (r) is com¬ 
pletely unconstrained by the Jeans analysis not only at 
the center (r = 0 ), but also at any other radius r within 
the stellar body of the dwarf. 

This paper is organized as follows. Section § [2] de¬ 
scribes the dSph model. Section § [3] discusses the obser¬ 
vational data (stellar line-of-sight velocities and stellar 
surface brightness profiles). Section § [4] gives a detailed 
description of my optimization technique. Section § [S] 
presents the main results of this study. The conclusions 
are presented in Section § [ 6 ] 

2 . MODEL 

Assuming a spherical symmetry, the Jeans equation for 
a two-component (DM -|- stars) system can be written as 

1 d(p^g^) ^ 4?7 cr^ ^ 

p* dr 1 + rj r dr 

Here r is the distance from the center, p*, Ui., and rj 
are the stellar density, radial velocity dispersion, and 
anisotropy, and $ is the total (DM -|- stars) gravita¬ 
tional pot ential. We introdu c ed the anisotropy param¬ 
eter p in iMashchenko et all (|2006all . It is defined as 
p = (ct^ — CT^)/(tT^ + 0 "?); and is related to the more con¬ 
ventional anisotropy parameter /3 through p = /3/(2 — /3) 
and /3 = 2p/(l -|- p). (Here at is the one-dimensional 
tangential velocity dispersion.) Unlike /?, the parameter 
p is conveniently symmetric: it is equal to — 1 , 0 , and 1 
for purely tangential, isotropic, and purely radial orbits, 
respectively. (The corresponding /3 values are —oo, 0, 
and 1 , respectively.) 

I assume that the stellar anisotropy p smoothly varies 
between the two asymptotic values, po at the center and 
Pi in the infinity: 


, . ^ poA + pi(r/ra)" 
’ A+(r/r,)2 


( 2 ) 


(|Baes &: van Hesell2007t iStrigari et all 1200811 . Here is 
the anisotropy transition radius, and A= (l-|-pi )/(!-!- 
Vo)- 

The DM halo is assumed to have the following density 
profile: 


/ 13-y - (3) 

{r/rs)^{l +r/rs)-^ ^ 

(Here rg and ps are the DM scaling radius and den- 
sity, respectively.) This is a generalization of the 
iNavarro, Frenk, fc Whit j (|1997l NFW) profile, with an 
arbitrary central logarithmic density slope 7 and the 
outer slope of 3. In the limit 7 = 1 we recover the theo¬ 
retical NFW profile; when 7 = 0, we obtain a flat-cored 
profile, which is a lmost identical to the observationally 
inferred iBurkertI (119951) profile. Stars are assumed to 
ha ve a generalized Plummer density profile (see eq. ([ 6 ]) 
in The DM contribution to the total gravitational 

potential $ is obtained by numerical integration of the 
corresponding Poisson equation; the stellar contribution 
has simple analytical solutions (different for different val¬ 
ues of the stellar outer density slope a). 

To compare the solution of the Jeans eq. 0 with ob¬ 
servations, one has to compute the stellar line-of-sight 


velocity dispersion, cxios, as a function of the projected 
radius, R: 


R 

(|Binnev fc Tremain^ll987L p. 208) Here E(i?) is the stel¬ 
lar surface density. 

Overall, the model has six free parameters: rg, Ps^ 7 
(characterize the DM distribution), and po, pi, (char¬ 
acterize the stellar anisotropy profile). 


1 - 




p*a;.r 


dr 


(U-p)rV (r2-i?2)i/2- 


( 4 ) 


3. OBSERVATIONAL DATA 
3.1. Surface brightness profiles 

For Jeans mass modeling it is essential to know stellar 
density profiles of dSphs as accurately as possible. These 
can be derived from surface brightness profiles if one as¬ 
sumes a spherical symmetry and a certain value for the 
stellar mass-to-light ratio. As Galactic dSphs are located 
within 300 kpc from the Sun and are fully resolved into 
stars, one has to employ star count techniques to derive 
their surface brightness profiles. 

Traditionally, a few simple analytical stellar sur¬ 
face density profiles have be en us e d for mass mod- 

eling of dSphs: empiric al iKind (11962 _ Plummer 

(jBinnev fc Tremaiii^ 119871). a nd Sersid ( 19631) . We 
showed in IMashchenko et al.l (l2006a|) that the surface 
density profile of the form 

S = So[l + (W]-^“-'^/^ (5) 

corresponding to the “generalized Plummer law”, 


p, = p,,o [1 + (r/b)'^] 


-ajl 


( 6 ) 


describes the surface brightness profile of Draco dwarf 
spheroidal extremely well. (Here R and r are the pro¬ 
jected and spatial distances from the dwarf’s center, Eq 
and are the central stellar surface and volume den¬ 
sities, b is the “core radius”, and a is an integer which 
should be ^ 4 for finite total stellar mass models.) The 
classical Plummer law is recovered when a = 4. 

For the five dSphs analyzed in this paper, I ob¬ 
tained the generalized Plummer profile parameters, b 
and a, by fittin g of eq. jsL t o the best avail¬ 
able star count dat a: iColeman et al.l (l2005bl their Fig¬ 
ure 15) for Fornax. iSmol cic et al .l ( 2007L their Figure 9) 
for Leo I, iGoleman et al.l (l2005al the i r Figu re 6) for 
Sculptor, and llrwlu fc Hatzidimitrioul (1199^ for Sex¬ 
tans. For Garina, there were significant discrepancies 
between the published star count profiles, so I car- 
ried out a joint fitt i ng for t hree different da t asets : 
lirwin fc Hatzidimitrioul (119951). iMaiewski et al.l (120001 
their Table 3). and iMunoz et ahl ( 20061 their Figure 4b). 
As the fitting was done in the linear space, all star 
count data points (even those which became formally 
negative after subtracting the assumed contribution from 
the background sources) were used. Figure [T] shows the 
star count data and the best fitting generalized Plummer 
models for the five dSphs. 

Table [T] presents the best x^ values (normalized by the 
number of degrees of freedom) for the five dSphs for four 
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log (R/b„.aj) 


Figure 1. Normalized radial surface brightness profiles. Point s with error bars show t he ob s erved star counts ; the lines s how the best 
fitting generalized Plummer models. For Carina, the data from ITrwin &: Hatzidimitrioul Ill995lh fMaiewski et al.l ll2QQg) . and lMuhoz et all 
(l20Qar) are shown as filled circles, empty circles, and diamonds, respectively. 


Table 1 

for different surface brightness models 


Galaxy 

Xgp 

Xf 

Xrk 

Xr 

Carina 

2.18 

2.18 

2.83 

2.92 

Fornax 

2.35 

9.81 

2.86 

2.55 

Leo I 

3.12 

47.8 

11.5 

6.66 

Sculptor 

3.50 

3.50 

4.86 

4.98 

Sextans 

0.87 

1.03 

0.84 

0.89 

Draco 

0.65 

1.69 

1.25 

1.10 

Average 

2.11 

11.0 

4.02 

3.18 


Note. — Xgp’ Xp? Xek’ xi normalized x^ values for 
the generalized Plummer, classical Plummer, empirical King, and 
Sersic profiles, respectively. 


simple analytical expressions - as long as a is an inte- 
ger) and appears to be the best fitting profile among 2-3 
parameter models. 

One complication is that all the published star count 
profiles were computed inside elliptical annuli (usually 
with the fixed ellipticity e and position angle), with 
the quoted radius measured along the major axis of the 
galaxy. For the purposes of this study, spherically aver¬ 
aged (in the plane of the sky) surface brightness profiles 
are required. For each galaxy, I computed the spher¬ 
ically averaged profile by solving numerically the fol¬ 
lowing integral for 100 values of the projected radius 
ii = 0.001... 30: 


different profiles (generalized Plummer, classical Plum¬ 
mer, empirical King, and S&sic); I al so included the fit¬ 
ting r esults for the Draco data from lOdenkirchen et ^ 
(120011 their samp l e S2) - the same data we used in 
iMashchenko et all (l2006all for mass modeling of this 
galaxy. As you can see in Table [TJ the generalized Plum¬ 
mer profile is the best one overall: for each galaxy, it 
produces the same or better fit than the other profiles, 
and is significantly better in terms of the average for 
all the six dwarfs. I warmly recommend this profile for 
all future dSph mass modeling work, as it is simple (sur¬ 
face and volume densities and gravitational potential are 


t/2 


51avr(^) — 


l-f 


i?2 


cos^ ip + [1 — e)2 sim p 


dp. 


(7) 

Here e is the ellipticity of the galaxy. The resultant, av¬ 
eraged profile is very close (but not identical) to the gen¬ 
eralized Plummer profile with the same exponent a and a 
smaller core radius b. I derived this spherically averaged 
value of b by fitting the above profile with the gen¬ 
eralized Plummer profile. The correction factor &/&maj 
is a function of both e and a, and ranges from 0.816 
(Sextans) to 0.895 (Leo I). (Here femaj is the core radius 
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obtained from fitting the generalized Plummer model to 
the original data, before averaging.) 

The results of the generalized Plummer profile fitting 
and other input galactic parameters for the five dSphs 
are presented in Table [2] As this study properly takes 
into account the self-gravity of stars, I had to assume 
certain values for the stellar mass-to-light ratios in V- 
band, (M/L)*. For each galaxy, I used the a verage of the 
Salpet er and composite (M/L)* estimates of iMateo et al.l 
(jl998l their Table 6 ). Knowing the distance to the 
galaxy, D, the F-band central surface brightness, Sq) 
the generalized Plummer model parameters a and b, and 
(M/L)* allowed me to compute the central stellar den¬ 
sity, (see Table [2]). The table also lists the estimated 
half-light radii (in arc minutes), i?hi, for spherically av¬ 
eraged models. 

3.2. Line-of-sight star velocities 

I used the stellar line-of-sight velocities catalogs of 
iWalker et al.l (|2009all for Carina, Fornax, Scu lptor, and 
Sextans, and the catalog of iMateo et al.l (|2008ll for Leo I. 
These are the largest homogeneous catalogs available for 
dwarf spheroidals, with the number of stars per galaxy 
ranging from ~ 300 (Leo I) to ~ 2500 (Fornax). I 
used the merged version of the datcQ. The median er¬ 
ror of the velocity measurements is ±2.1 km s“^. In 
all the dwarfs (except for Leo I) some contamination 
from foregro und Galactic s t ars is e xpected to be present 
in the data. IWalker et al.l (|2009bD quantified this con¬ 
tamination by computing individual “membership prob¬ 
abilities” P for each star. For this study, I only used 
stars with P > 0.95 for Carina, Sculptor, and Sextans, 
and stars with P > 0.75 for Fornax. (Despite the less 
conservative selection criterion for Fornax, the resultant 
line-of-sight dispersion profile looks very smooth and un¬ 
contaminated, and the normalized for ttie best fitting 
model is ^ 1 - see Figure O) Due to the Leo Ps large 
heliocentric velocity and distance from the Sun, the data 
for this dwarf should have essentially zero contamina¬ 
tion in the r ange of heliocentric velocities from 260 to 
315 km s“^ (| Walker et al.ll200'^ their Fig. 1). Accord¬ 
ingly, I used these simple velocity cuts to select the Leo I 
star candidates. 

Carina was the only galaxy in our sample for which 
the best fitting model deviated more than two sigmas 
from the data. This could be due to unaccounted fore¬ 
ground contamination (which would not be surprising, 
given that it is the closest to the Galactic plane dSph 
and has a very low surface brightness). For this galaxy 
only, I removed the six stars (after removing the stars 
with P < 0.95), for which the deviation from the best 
fitting model was more than three sigmas, and then re¬ 
peated the full mass modeling for the reduced star list. 
(Here one sigma corresponds to the total, model ± obser¬ 
vational uncertainty, velocity dispersion.) The removed 
stars ha ve the following designations in the IWalker et al.l 
(I2nn9all catalog: Car-0056, 0200, 0465, 0547, 0680, 1892. 
This operation boosted the absolute probability of the 
best fitting model for Carina from 0.04 to 0.26 (when 
A^bin = 7; see ? 14.31) . It is important to note that all the 
six stars are located in the outer half of the radial bins, 
so this operation should have had very little impact on 

^ http://www.astro.Isa.umich.edu/$\sim$nmiat6o/research.html 


the analysis of the central parts of the galaxy - which is 
the main focus of the current study. 

4. CHOOSING GOOD MODELS 
4.1. fitting 

Two popular statistical methods used in mass mod¬ 
eling of dSphs are fitting and maximum likelihood 
analysis. Both methods have advantages and disadvan¬ 
tages. The x^ method complements naturally the Jeans 
analysis, as both operate only with the second moment 
(dispersion) of the stellar velocity distribution, whereas 
for the maximum likelihood approach one has to assume 
a certain (usually Gaussian) form for the line-of-sight 
velocity probability distribution function. This factor 
makes the x^ fitting less sensitive to outliers, and hence 
a much better approach when one is interested in ab¬ 
solute probabilities of different models. The maximum 
likelihood method, on the other hand, is more appropri¬ 
ate for post-Jeans mass modeling, where the full model 
line-of-sight velocity PDF is computed and compared to 
the data. 

The main disadvantages of the x^ technique are that 
1 ) one has to try a few different values for the number of 
the radial bins, A^bin ~ but that only means longer com¬ 
puting time, and 2) in its simplest implementation, the 
random locations of individual stars within each radial 
bin is an additional source of variance in the analysis. To 
address the latter issue, for both the observational data 
and model I estimate the following quantity (for each 
radial bin): 



i=l 


Here ai is the line-of-sight velocity dispersion at the spe¬ 
cific location (projected distance from the dwarf’s cen¬ 
ter) of the f-th star, and n is the number of stars in the 
current radial bin. The quantity S will converge to the 
square of the local line-of-sight velocity dispersion when 
the radial extent of the bin shrinks to zero and the num¬ 
ber of stars (n) grows to infinity. The important point 
is that such a convergence is not required when all one 
wants to do is to measure the statistical deviation of the 
model from the data. All is needed is to estimate the 
mean (for the data and the model) and standard devia¬ 
tion (only for the data) of the quantity S for each radial 
bin, and then compute the x^ deviation for S from all 
the radial bins. As the S measurements for both ob¬ 
servational data and model are carried out for the same 
projected distances from the center of the dwarf (corre¬ 
sponding to the distances of actual stars), the variance 
due to random locations of stars within radial bins is 
completely removed. 

Of course, we do not know from the observations the 
values of ai at the position of each star. Instead, we know 
the individual line-of-sight velocities (corrected for the 
systemic velocity of the galaxy), F, and the associated 
measurement uncertainties, Cmes.i- One can use these 
known quantities to write down the following unbiased 
estimator of S (without any assumptions on the shape of 
the velocity PDF): 
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Table 2 

Input galactic parameters 


Galaxy 

D 

kpc 

So 

mag arcsec ^ 

e 

(M/L)* 

M 0 /L 0 

V"sys 

km 

a 

^maj 

arcmin 

b 

arcmin 

p*,0 

Mq pc“^ 

^hl 

arcmin 

Garina 

101 

25.5 

0.33 

0.845 

222.9 

5 

10.89 

9.00 

0.0055 

TFT 

Fornax 

138 

23.4 

0.31 

0.935 

55.2 

7 

27.92 

23.61 

0.0147 

7.67 

Leo I 

255 

22.6 

0.21 

0.74 

282.9 

13 

10.21 

9.14 

0.0490 

1.92 

Sculptor 

79 

23.7 

0.32 

1.23 

111.4 

5 

9.95 

8.28 

0.0584 

3.66 

Sextans 

86 

26.2 

0.35 

1.15 

224.3 

4 

22.18 

18.09 

0.00195 

10.45 


Note. — Here D is the distance from the Sun, So is the central surface brightness in y-band, e is the ellipticity, (M/L)* is the stellar 
mass-to-light ratio in V-band, \4ys is the systemic heliocentric velocity, a and bmaj th® generalized Plummer parameters, b is the 
spherically averaged generalized Plummer par ameter , p*_n is the central stellar density, and /?hl is the half-light radius for a spherically 
averag ed galaxy. D, Sq, and e were t a ken fr om lMat^ l|199gl) for all the galaxies except for Leo I; for Leo 11 used the data from IMateo et ^ 
1120081) . (M/L)* is from IMateo et ahl 1119981 . their Table 6). Vsys is from lWalker et all ll2009bl) . The rest of the parameters were derived in 
this paper. 


1 ^ 

(9) 

This expression is valid even for the cases when the line- 
of-sight velocity dispersion changes dramatically within 
a single radial bin. The standard deviation of S can be 
estimated from the following expression, 



IL 2 

2 E 


( 10 ) 


which does require the velocity dispersion to be approx¬ 
imately constant within a single bin. Numerical experi¬ 
ments showed that eq. m noticeably deviates from the 
true answer only when the velocity dispersion changes 
significantly (by a factor of two or more) within a sin¬ 
gle radial bin, which is never the case in my best fitting 
models. 

Eqs. (l 8 UIO)l are written for a single radial bin. It is now 
straightforward to write down the expression for the 
deviation of the model from the data for all the Vbm 
radial bins: 


-^bin 

= E 

i=i 



( 11 ) 


Here Sj , Sj , and asj are computed for the j-th radial bin 
using eqs. and OT, respectively, and k = Nhin — 

A is the number of the statistical degrees of freedom, 
where A the effective dimensionality of the free parameter 
space (see ? 14.21) . 

The only drawback of the above procedure is that 
it can be substantially more computationally expensive 
than the standard analysis, as one has to compute 
the model line-of-sight velocity dispersion (numerical in¬ 
tegration using eq. @) for each model, each star, and 
each value of iVbin- 


4.2. Artificial data tests 

One can already use eq. m from the previous section 
to compute the relative goodness of fit for different mod¬ 
els. To go one step further and convert the x^ values to 
the absolute probabilities of the models. 


n = p(if). (12) 

one has to estimate the effective number of dimensions 
of the free parameter space, A. (Here P{x,y) is the reg¬ 
ularized Gamma functior0-) 

To estimate A for my six-parameter model described 
in § [21 I generated a dataset of 2499 fake stars located 
at the same projected distances from the dwarf’s center 
as the real 2499 stars in Fornax with known line-of-sight 
velocities (see §IS21). I used my best fitting Fornax model 
(for the iVbin = 7 case) and actual individual star veloc¬ 
ity observational uncertainties, CTmesj to generate random 
values of the “observed” line-of-sight velocities for the 
artificial stars. (I assumed that both the model line-of- 
sight velocities and the observational velocity uncertain¬ 
ties have Gaussian distributions.) 

Next, I split my artificial dataset into 21 equal size sub¬ 
sets (each one consisting of 119 randomly selected fake 
stars). In a sense, I generated 21 random realizations of 
the Fornax galaxy. Then I carried out full mass modeling 
(as described in § 14..'ll with A^bin = 7) and derived the 
xl value (eq. (ITT|) 1 for the best fitting model, separately 
for each sub-set. Assuming that the resultant 21 random 
numbers are drawn from the same Xk distribution, their 
mean should be an unbiased estimator for k. In this nu¬ 
merical experiment, I derived k ~ 3.73, corresponding to 
the effective number of dimensions in my free parameter 
space A = A^bin — A: ~ 3.27. I used this number in the 
rest of the paper. This number is significantly (almost 
by a factor of two) lower than the formal number of the 
model free parameters (six) reflecting significant degen¬ 
eracies (from the Jeans analysis point of view) present 
in the model. As one will see in § |5l the most degen¬ 
erate is the parameter 7 (the logarithmic slope of the 
DM density). Post-Jeans mass modeling (dealing with 
the full line-of-sight velocity PDF, and not just with the 
dispersion) is expected to be less degenerate. 

It is clear that A will not be exactly the same for differ¬ 
ent data and for different number of radial bins (because 
the radial extent of the binned data depends on N\^in), 
but as long as the number of bins is relatively large (so 
that k = A^bin — A 1), the computed absolute proba¬ 
bilities should not be very sensitive to the errors in A. To 

^ http: //en. Wikipedia. org/wiki/Chi-square_distribution 
























Mass Models of Dwarf Spheroidal Galaxies 


7 


be more quantitative, I computed numerically the rela¬ 
tive error for the absolute probabilities If resulting from 
the A uncertainty, (AA/II) dll/dfc, for k = 3.73 and two 
different cases: = k (corresponding to 11 ~ 0.4, as in 

my best fitting models), and xi ~ 2.5fc (corresponding 
to n ~ 0.04, my lower cutoff for “good models”). Purely 
statistical (Poissonian) errors in estimating k (and hence 
A) from the limited sample of the size Ns = 21 is of 
the order of {2k/NsY^'^ — 0.6. Assuming a larger uncer¬ 
tainty (to account for a variance in data) A A = 1, the 
resulting relative errors in 11 are only 0.05 and 0.35, re¬ 
spectively, which I consider to be acceptable. If Abin > 7, 
as is the case with most of my galaxies, the relative errors 
in n should be significantly smaller. 

4.3. Numerical algorithm 

In this section I outline my overall mass modeling pro¬ 
cedure. 

I thoroughly explored (“brute force optimization”) the 
four out of the six free model parameters: the central 
logarithmic DM density slope 7, and the three stellar 
anisotropy parameters, 770, ryi, and (see §[ 21 )- For 7, I 
used the discrete values ( 0 , 0 . 1 , 0 . 2 , ..., 1 ), which cov¬ 
ered the whole plausible range - from the observations- 
suggested flat core (7 = 0) to the CDM prediction, 
7 ^ 1 . In these low stellar density objects the presence 
of baryons is not expected to make the DM distribution 
“cuspier” via the adiabatic contraction of DM. A pres¬ 
ence of a sufficiently massive central black hole could 
have steepened the DM cusp, but there is no evidence 
for super-stellar mass black holes in dSphs (iLora et al.l 

IMi). 

Both rjo and rji had the following discrete values: (—1, 
—0.9, —0.8, ..., 1), which covered the full range of 
anisotropies - from purely tangential (77 = — 1 ) to purely 
radial (77 = 1). In addition, for the central anisotropy rjo I 
only explor ed the range a l lowed by the central anisotropy 
theorem of I An fc EvansI (120061 1. 770 < ( 4/7 — 1)“^. The 
anisotropy radius had the following values (in the units 
of the stellar core radius b): (0.1, 0.147, 0.215, 0.316, 
0.464, 0.681, 1.47). They are equidistant in the loga¬ 
rithmic space, and span the full range of radial distances 
resolved by the data. Overall, I explicitly tested 22,344 
different models with the above discrete values of the four 
model parameters. 

The remaining two model parameters, DM halo scaling 
radius and density ps , are very well constrained by the 
data when the rest of the parameters are fixed, which al¬ 
lowed me to use implicit optimization technique for these 
two parameters. This dramatically reduced the number 
of models to be tested and hence the total computational 
time, as rg and especially ps span a very large range of 
values (four orders of magnitude for ps, see Table [3]). 
The only prior was for rg to stay within the interval 0.05 
... 5 (in 6 units; for Sculptor, I used a different range 0.1 
... 15, as this galaxy has a very compact core). There 
was no prior for pg. 

I used the one-dimensional Brent’s method 
(jPress et al.lll992ll separately for rg and pg optimization 
(finding global minimum in pg optimization is the 
inner loop, and rg optimization is the outer loop). The 
method showed a good numerical convergence, requir¬ 
ing only of order of a hundred models with different 


combinations of (rg, ps) to be tested (when the other 
four explicit parameters were fixed). Including both 
explicitly and implicitly explored parameters, I had to 
compute the deviation from the data for ~ 2 x 10 ® 
different models - per galaxy and per Abin value. 

For a given galaxy and a chosen combination of the 
six model parameters and the value of Abin, I solved the 
Jeans eq. o numerically, and then solved numerically 
the integral in eq. ( 0 ]) to derive ctios at all the projected 
radii R corresponding to the locations of the dwarf’s stars 
with known line-of-sight velocity. Then I used eqs. da- 
El) to compute the deviation of the model from the 
data, and eq. ( 11 ^) to assign absolute probabilities to the 
models. 

Absolute probabilities of the best fitting models, Ilmax, 
range from (1.23 for Sculptor to 0.60 for Sextans (Table|31), 
corresponding to ~ 0.5 —1.2 sigma deviations. I used one 
global lower cut-off value of Blow = 0.0455 (a two sigma 
deviation) to select good models for all the galaxies. 

I used the following approach to determine the optimal 
value of the number of radial bins, Abin, for each galaxy. 
In the limit of large Abin (small number of stars in one 

bin, n), as oc N^^^ (see ea. ifTUll l and the x^ analysis 
becomes dominated by the Poissonian noise and hence 
“fuzzier” (less discriminatory). In the opposite limit, 
Abin — 0 , the signal to noise ratio improves, but the x^ 
analysis becomes again less accurate because with fewer 
radial bins we start loosing radial resolution and may 
overlook some sharp ctios features near the dwarf’s cen¬ 
ter or in the galactic outskirts. In addition, when Abin 
approaches A (the effective dimensionality of the free pa¬ 
rameter space), the errors in measuring A become more 
important, which makes the absolute model probabilities 
less reliable (see S I4.2|1 . Hence it appears there should be 
a certain value of Abin when the x^ analysis is the most 
accurate (or discriminatory). I found this optimal value 
separately for each galaxy by running the full set of mod¬ 
els for Abin = 7 (minimum acceptable number of bins for 
this model, see § [ 121 ), 15, and 30, and then choosing the 
value which produced the smallest fraction of good mod¬ 
els (/ 2 cr in Tabled]). For Fornax, Sculptor, and Sextans 
the optimal value of Abin was found to be 15; for Carina 
and Leo I it is equal to 7 (Table |3|). 

5. RESULTS 

I applied the mass modeling procedure described in the 
previous sections to all the hve dwarfs. The full results 
are available online. As good models found in the cur¬ 
rent analysis occupy only a small fraction of the total 
six-dimensional free parameter space, the data presented 
here can be used for post-Jeans mass modeling projects 
(utilizing the full PDF of the observed stellar line-of-sight 
velocities; the subject of Paper II) to dramatically reduce 
the required computational time. 

Figure [2] shows the radial uios profiles for the two best¬ 
fitting models (with 7 = 0 , solid lines, and 7 = 1 , dotted 
lines) for each galaxy. As one can see, the overall quality 
of the model hts is excellent, with the x^ values (normal¬ 
ized by the number of degrees of freedom k = Abin — A) 
being close to a unity. It is remarkable that these two 
dramatically different values of 7 (0 and 1 ) make very lit¬ 
tle difference in the quality of the model fit. This is most 
obvious for the galaxy with the highest quality data - 












Mashchenko 



Figure 2. Radial line-of-sight stellar velocity dispersion profiles. Solid circles with error bars show the binned observational data. The 
solid and dotted lines show the best models with 7 = 0 and 1 , respectively. 


Fornax - where the two ctios profiles are essentially iden¬ 
tical way beyond the radial range covered by the data. 
For the galaxies with the best 7 = 0 and 7=1 profiles 
apparently diverging near the center (Sculptor and Sex¬ 
tans), one could naively think that with a larger iVbin 
(and hence with the binned data covering a wider range 
of radial distances) the two models would become more 
statistically distinct. But this is not the case: for all the 
values of iVbin I tried (7, 15, and 30), and for all the five 
galaxies, the difference between the best fitting models 
with different values of 7 is statistically insignificant, to 
a similar degree. 

Tables [3] and |4] summarize the main results of this 
study. In particular, I list there the ranges spanned by 
good (< 2a) models separately for the free model param¬ 
eters rjo, Vs, and ps- As one can see, these parameters are 
very poorly constrained by the Jeans anal ysis. The full 
allow ed range for ryo is —1 ■ • ■ 0.3 (due to the I An fc EvansI 
120001 central anisotropy theorem), and only three galax¬ 
ies (Fornax, Leo I, and Sculptor) have a (very weak) 
constraint on this parameter, in the sense that the pure 
central tangential anisotropy is excluded. Similarly, only 
one galaxy (Sculptor) has any constraints on the DM 
halo scaling radius: Vs > 0.5 (in b units). For the rest of 
the dwarfs, good models span the full explored range of 
Cg, from 0.05 to 5. Parameter ps is the only one in the 
current analysis which is allowed to take any values (no 
priors), so it does appear to be constrained in a mean¬ 


ingful way. One has to remember though that there is 
a strong degeneracy between ps and Vg. E.g., if were 
allowed to go to even smaller values (< 0.05), the range 
for Ps would expand to even larger values. I do not list 
global constraints for the three remaining model param¬ 
eters, 7 , 71 , and Tq, as good models span the full ranges 
for these parameters: 0 ... 1, —1... 1, and 0.1... 1.47, re¬ 
spectively. The conclusion here is that in mass modeling 
of dwarf spheroidal galaxies the Jeans analysis is not ca¬ 
pable of placing meaningful constraints on global model 
parameters (for models which are flexible enough, as in 
the current study). 

The analysis does constrain the minimum central DM 
density, po.min, and the corresponding mass to light ratio, 
(M/L)o,min (see Table [3]). Sculptor and Leo I have the 
strongest constraints on po.min (^ 0.09 Mq pc“^ and 
> 0.08 Mq pc“^, respectively, at the two-sigma level). 
Carina, on the other hand, has the strongest (M/L)o,min 
constraint: > 8 Mq/Lq. 

More useful constraints can be derived for radius- 
dependent (computed at different distances from the 
dwarf’s center) quantities. Figure [3] plots the range of 
the averaged DM density pavr spanned by good mod¬ 
els as a function of the enclosed radius r, separately for 
each galaxy. The averaged stellar density profiles (dash- 
dotted lines) are shown for comparison. There has been a 
debate in the literature (IMateo et al.lll993l : IStrigari et al.l 
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Table 3 

Global constraints 


Name 

N, 


f2a 

Hmax 

Vo 


Ts 


ps 

PO.min 
Mq pc-3 

(M/L)o,min 

-(^0 A© 

Carina 

740 

7 

0.680 

0.262 

-1... 

0.3 

0.05., 

, .5 

1.7... 2.5 X 10"" 

0.048 

8.2 

Fornax 

2499 

15 

0.253 

0.468 

-0.8.. 

.0.3 

0.05., 

, .5 

0.25... 2.1 X 10^ 

0.028 

2.7 

Leo I 

328 

7 

0.210 

0.302 

-0.8.. 

.0.3 

0.05., 

, .5 

0.13... 600 

0.080 

1.95 

Sculptor 

1349 

15 

0.230 

0.228 

-0.9.. 

.0.3 

0.5... 

15 

0.16.. .63 

0.087 

3.1 

Sextans 

397 

15 

0.906 

0.598 

-1... 

0.3 

0.05., 

, .5 

0.90... 2.9 X lO'^ 

0.007 

5.3 


Note. — Here A/"* is the number of stars with a known line-of-sight velocity used for the analysis; A^bin is the number of radial bins; 
/ 2 cr is the fraction of the “good models” (with the absolute probability H > 0.0455); Hmax is the highest absolute probability for a model. 
I also list the two-sigma ranges spanned by good models for the global model parameters tjq (central anisotropy), Vs and ps (DM scaling 
radius and density). Finally, the two-sigma constraints on the minimum central DM density (po,inin) and the corresponding mass to light 
ratio ((M/L)o,min) are given. 


Table 4 

Radius-dependent constraints 


Name 

ri 

pc 

Ml 

10’’Mq 

(M/L)l 

r 2 

pc 

Pdm 

Mq pc“® 

Mq/Lq 

73 

pc 

V 

M 300 

IO'^Mq 

Meoo 

WMq 

Carina 

409 

0.74... 

1.14 

26. 

. .39 

187 

0.032. 

. . 0.074 

14. 

. .32 


-1.. 

.0.3 

0.35.. 

.0.78 

0.85. 

.2.3 

Fornax 

924 

6.1... 

8.0 

5.5. 

. .6.9 

453 

0.021. 

. . 0.040 

3.7. 

. .6.2 

232 

-0.88. 

. 0.05 

0.27. 

.3.4 

1.90. 

.5.1 

Leo I 

389 

1.55. .. 

1.94 

3.4. 

. .4.1 

204 

0.067. 

. .0.115 

2.5. 

. .3.8 

72 

-0.67. 

..0.26 

0.76.. 

. 1.43 

1.94. 

.5.7 

Sculptor 

435 

2.17... 

3.01 

22. 

. .30 

256 

0.067. 

. .0.107 

20. 

. .31 

83 

-0.46. 

..0.19 

0.79.. 

. 1.63 

00 

16 

Sextans 

1034 

1.4... 

3.2 

19. 

. .42 

451 

0.004. 

. .0.013 

11. 

. .32 


-1.. 

.0.3 

0.06. 

. 1.1 

0.37. 

. 1.7 


Note. — Here ri is the radius where the enclosed DM mass has the tightest constraints. At this radius, I list the two-sigma constraints 
on the enclosed mass (Mi) and the averaged light to mass ratio ((M/L)i). Similarly, r 2 is the radius where the local DM density has 
the tightest constraints, with the corresponding two-sigma constraints on the local DM density (pdm) and the local mass to light ratio 
((M/L) 2 ). Finally, is the radius where the local stellar velocity anisotropy, 77 , has the tightest constraints. In addition, I show the 
two-sigma constraints for the enclosed mass at the radii 300 pc (M 300 ) and 600 pc (Mgoo)- 


120071 [2008tlWalker et al.l[2009cl : lWolf et al.l[2009D regard¬ 
ing the value of the radius at which the Jeans analysis 
provides the tightest constraints on the enclosed mass 
(or equivalently the averaged density). As one can see in 
Figure [21 such radius indeed exists (vertical dotted lines) 
and is very well defined. In TablejUl list these radii (ri), 
the range of the enclosed mass spanned by good models 
at those radii (Mi), and also the corresponding range 
of the mass-to-light ratio ((M/L)i), for all the dwarfs. 
For comparison, in Table |2] I also list the good model 
ranges of the enclosed mass at 300 and 600 pc - M 300 
and Mgoo, respectively. As one can see, ri ranges widely, 
from 390 pc for Leo I to more than 1000 pc for Sextans, 
and at neither of the previously suggested radii (300 and 
600 pc) the enclosed mass is well constrained by the Jeans 
analysis for all the dwarfs. To be more quantitative, the 
mean-squared two-sigma error for all the five dwarfs is 
±0.40 dex, ±0.24 dex, and ±0.10 dex (±25%) for the 
enclosed mass measured within 300 pc, 600 pc, and ri 
(unique for each galaxy), respectively. The accurately 
determined values of the enclosed mass presented here 
can be very useful for identifying the dSph counterparts 
in cosmological simulations, and is probably the most 
important result of t he current study. 

The conjecture of iWalker et al.l (l2009cll that the en¬ 
closed mass is well constrained at the half-light radius, 
i?hh is not corroborated by the current study. In fact, en¬ 
closed masses measured at this radius result in the worst 
mean-squared two-s igma error, ±0.58 dex. 

The conclusion of lWolf et al.l (|2009ll (whose analysis is 
hard to interpret, as DM and stars are coupled in their 
model) that the enclosed mass has the best constraint 


at the radius where 7 * = —dlogp*/dlogr = 3 is also 
not borne out in the current study: for the five dwarfs 
I obtained 7 * = 3.54 ± 0.38 at my best radii (ri), with 
the full range from 3.22 to 4.20. My sample of galaxies 
is too small to try to find a good estimator for ri, but I 
believe it should be a function of both the stellar density 
profile ( 7 *(r)) and the quality and quantity of the stellar 
line-of-sight velocity measurements. This issue needs to 
be explored further once a larger sample of dwarfs is 
subjected to a similar analysis. 

To estimate how sensitive the above analysis is to par¬ 
ticular values of fVbin, I repeated it for all the three val¬ 
ues of Vbin: 7, 15, and 30. For the less optimal values 
of Vbin, the radii where the enclosed mass is measured 
to the highest possible accuracy (r,) does not deviate 
from the value obtained with the optimal Vbin by more 
than 7% (for Sextans; for Fornax it is essentially zero). 
The averaged DM density measured at ri is also very 
close for all the values of A^bin (less than one sigma de¬ 
viation), with one notable exception: for Leo I, moving 
from iVbin = 7 (the optimal value) to Abin = 30 re¬ 
sults in pavr(?'i) becoming smaller by 0.18 dex, or 3.8 ct. 
(This may be related to the fact that Leo I has the small¬ 
est number of stars with a known line-of-sight velocity: 
A* = 328.) Overall, this analysis seems to be largely in¬ 
sensitive to how well the optimal value of Abin is chosen. 

Another proposition discussed in the literature is that 
there is a radius (300 pc) at which the enclosed mass 
is the same for all dwarf spheroidals - around 10 ^ Mq 
([M ateo et al.lll993t iStrigari et al.l 12008( 1. This is barely 
(at 2(7 level) consistent with the data in the current anal¬ 
ysis, as there is a tension between Carina on one side, and 
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Figure 3. Averaged DM density, pavr, as a function of radius, r, in stellar units (p*,o and 6, respectively). Two solid lines outline the range 
of the models consistent with the observational data at better than 2cr level. Dashed lines correspond to formally best models. Vertical 
dotted lines mark the radius (ri in Table[3ll where the pavr uncertainty is the smallest. Dash-dotted lines show the averaged density profiles 
for stars. 


Leo I and Sculptor - on the other side (see Table E]). The 
galaxies may have a similar mass at this radius, but it 
is unlikely to be identical. If one considers even smaller 
radii (say, 200 pc), then formally there will be a com¬ 
mon enclosed mass for all the five dwarfs consistent with 
the data at better than two sigma level, but only be¬ 
cause the uncertainties in measuring the enclosed mass 
increase dramatically at smaller radii (see Figure E]). 

Another interesting radius-dependent quantity, local 
DM density, is more model dependent, and is not as 
well constrained by the Jeans analysis as the enclosed 
mass. Figure m shows the ranges for puM{r) spanned by 
good models. I also plot there the stellar density pro¬ 
files, p*(r), as dash-dotted lines. As with the averaged 
density case, each galaxy has a well defined radius, r 2 , 
(vertical dotted lines) where the local DM density is best 
constrained by the data. Table E] lists the values of r 2 
and the corresponding ranges of pdm(?' 2 ) and the local 
mass-to-light ratio {M/L )2 for each galaxy. r 2 tends to 
be a factor of two smaller than r,. The uncertainty in 
Pdm(?' 2 ) ranges from 0.20 dex for Sculptor to 0.5 dex 
for Sextans, which is substantially worse than the con¬ 
straints on the enclosed mass. Still, knowing the local 
DM density at a certain radius (which happens to be 
quite small) to a reasonably good accuracy can be very 
useful - for example, for studies aimed at direct search of 


DM in dwarf spheroidals via its annihilation signal. Both 
Figures Eland El clearly demonstrate that these dSphs are 
DM dominated at all radii. 

The last radius-dependent quantity which has useful 
(albeit very weak) constraints is the local stellar velocity 
anisotropy, r]{r). Figure El demonstrates that better than 
two sigma models span a wide range of rj at each radius 
r, which is only slightly narrower than the full allowed 
range (the space between the two dotted lines). The 
best radius, ra, is not as well defined as with /5avr(?’) and 
Pdm(»")- I list the values of and the corresponding 
good model ranges for r]{r 3 ) in Table El Sculptor has 
the tightest constraints on 77 (^ 3 ); Carina and Sextans are 
essentially unconstrained. 

Given how poor the constraints on r]{r) are, very lit¬ 
tle can be learned from the analysis of Figure El One 
interesting conclusion is that at any given radius the 
data are consistent with stellar velocities being isotropic 
{r]{r) = 0). One cannot conclude from this radius- 
dependent analysis that the global stellar isotropy is also 
consistent with the data. To address the latter issue, one 
has to analyze the full results of the current study (avail¬ 
able online). Specifically, one has to test if all the globally 
isotropic models (with rjo = rji = 0 ) are good ones (de¬ 
viate less than 20 ^ from the data). I carried out such an 
analysis, with the conclusion being that indeed the data 


























Mass Models of Dwarf Spheroidal Galaxies 


11 



Figure 4. Local DM density, p, as a function of radius, r, in stellar units (p*,o and b, respectively). Two solid lines outline the range 
of the models consistent with the observational data at better than 2a level. Dashed lines correspond to formally best models. Vertical 
dotted lines mark the radius (r 2 in Table[3ll where the p uncertainty is the smallest. Dash-dotted lines show the density profiles for stars. 


is consistent (at better than 2a level) with the global 
stellar velocity isotropy hypothesis f or all the dwarfs. 

It has already bee n demonstrated (|Strigari et al.ll20d8l : 
I Walker et al.ll200^ that in dSphs the central logarith¬ 
mic DM density slope, 7 , is not constrained by the 
Jeans analysis (when no stellar proper motion data are 
available). My analysis generalizes this result for all 
radii r: the local logarithmic DM density slope, 7 (r) = 
—dlogpDM/dlogr, is found to be unconstrained by the 
Jeans analysis. 

6 . CONCLUSIONS 

I carried out the classical Jeans mass modeling of the 
five dwarf spheroidal galaxies with the highest quality ob¬ 
servational data - Carina, Fornax, Leo I, Sculptor, and 
Sextans. My primary goal was to push the analysis to 
its limits, by developing a flexible enough galactic model 
(with variable stellar velocity anisotropy), refining the 
traditional fitting algorithm, and investing a signifi¬ 
cant amount (3 x 10® cpu hours) of supercomputing time 
to carry out an exhaustive model parameter optimiza¬ 
tion. The main results of this study are as follows. 

• My galactic model with the six free parameters (ts , 
Ps, 7 , rjo, 71 , and r^) gives a good description of the 
observational data: the normalized for the dif¬ 
ferences between the model and observed velocity 


dispersion profiles are close to a unity, which is con¬ 
sistent with the deviations being purely due to the 
observational errors. 

• I show that the Jeans mass modeling approach 
(even with high quality data and modeling) can¬ 
not place meaningful constraints on most galactic 
parameters, including the central DM density loga¬ 
rithmic slope ( 7 ), DM halo scaling radius (r^) and 
density (ps), and central stellar velocity anisotropy 
{rjo). Moreover, I show that the local DM density 
logarithmic slope, 7 (r) = —dlogpDivi/dlogr, is un¬ 
constrained at all the radii probed by the data. As 
a consequence, one must resort to much more com¬ 
putationally expensive post-Jeans mass modeling 
techniques (which deal with the full probability dis¬ 
tribution function for stellar line-of-sight velocities, 
and not just the dispersion) to be able to solve the 
“cusp - core” issue for dSphs, or to identify the¬ 
oretical mechanisms responsible for flattening DM 
cusps in dwarf galaxies (by constraining the stellar 
velocity anisotropy profiles). 

• My most important finding is that there is a cer¬ 
tain radius (ri; different for each galaxy) where 
Jeans mass modeling provides the tightest con¬ 
straints on the enclosed mass in dSphs. For Ca- 
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Figure 5. Stellar anisotropy, 77 , as a function of radius, r (in units of b). Two solid lines outline the range of the models consistent with the 
observational data at better than 2cr level. Two dotted lines show the total allowed range of models. Dashed lines correspond to formally 
best models. Vertical dotted lines mark the radius (rs in Table [Sj where the 77 uncertainty is the smallest. 


rina, Fornax, Leo I, Sculptor, and Sextans my two- 
sigma constraints for the enclosed DM mass are 
0 . 74 ... 1.14 X 10^ Mq (at 409 pc), 6.1...8.0 x 
10^ M 0 (at 924 pc), 1.55... 1.94 x 10^ Mq (at 
389 pc), 2.17...3.01 x lO’’ Mq (at 435 pc), and 
1.4...3.2 X lO”^ Mq (at 1034 pc), respectively. 
These constraints are much tighter than the con¬ 
straints at the previously suggested “good” radii 
- 300 pc, 600 pc, and the half-light radius. The 
tight constraints can be very valuable for placing 
the dSphs in the proper cosmological context. 

• I also derive useful constraints on the local DM den¬ 
sity. Similarly to the enclosed mass, each galaxy 
has a certain radius (r 2 ) where my Jeans anal¬ 
ysis places the tightest constraints on this quan¬ 
tity. For example, 1 show that in Sculptor poM = 
0.067... 0.107 Mq pc“^ (two-sigma interval) at 
256 pc from its center. Also, the analysis pro¬ 
duced useful minimum central DM density (po.min) 
constraints (the largest one is for Sculptor: po ^ 
0.09 Mq pc“^). These constraints can be used in 
projects aimed at detecting DM in dSphs via its 
annihilation signal. 

• I show that stellar anisotropy profiles are very 
poorly constrained in Jeans mass modeling. The 


only useful results here are that the pure central 
tangential anisotropy is excluded at better than 
two-sigma level for three out of the five dSphs, and 
that the data are consistent with the global stellar 
velocity isotropy for all the five dwarfs. 

• A significant advantage of the exhaustive search 
through the multi-dimensional model parameter 
space employed in this study, with all the interme¬ 
diate results stored and available online, is that the 
results can be used for many other projects. Most 
importantly, this analysis can be used as the first 
preliminary step (which eliminates a vast majority 
of models which are incompatible with the data) in 
post-Jeans mass modeling projects (which model 
the full PDF for the stellar line-of-sight velocities) 
- the subject of my Paper II in this series. The 
latter approach should be able to overcome most 
of the degeneracies present in Jeans mass modeling 
(as exposed in the current work), hopefully settling 
once and for all whether dSphs have flat DM cores. 


The simulations reported in this paper were car¬ 
ried out on facilities of the Shared Hierarchical 
Academic Research Computing Network (SHARC- 
NET:www.sharcnet.ca). This research was supported in 
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